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SECTION  I 


DISCUSSION  OF  THE  VARIOUS  TRANSONIC  SMALL  PERTURBATION 
EQUATIONS  DESCRIBING  THREE-DIMENSIONAL 
WING- BODY  FLOWS 


1 . INTRODUCTION 

One  of  the  most  difficult  problems  yet  to  be  solved 
to  a certain  degree  of  completion  is  the  three-dimensional 
transonic  interference  problem.  At  transonic  Mach  numbers, 
three-dimensional  effects  due  to  the  presence  of  a fuselage 
or  wing  tip  can  influence  the  flow  field  for  large 
distances  in  the  spanwise  direction.  For  low  and  moderate 
aspect-ratio  wings  combined  with  a body,  the  flow  fields 
are  rarely  similar  to  two-dimensional  ones  anywhere  on  the 
wing  except,  perhaps,  on  certain  localized  regions. 
Therefore,  the  treatment  of  the  three-dimensional  transonic 
flow  problems  involves  more  than  a trivial  extension  of 
existing  two-dimensional  methods,  which  are  fairly  well 
developed  by  now. 

Some  of  the  major  complications  are: 

(1)  The  circulation  changes  with  the  spanwise 
coordinate  and  has  to  be  determined  as  part  of  the 
solution. 

(2)  The  swept  and  tapered  planform  shapes  complicate 
application  of  wing  boundary  conditions  using 
finite-difference  methods. 

(3)  The  shocks,  which  are  nearly  normal  to  the 
free-stream  direction  in  two  dimensions,  can  be 

"oblique  in  the  lateral  direction  for  finite 
swept  wings,  thereby  making  them  difficult  to 
capture  sharply  because  of  poor  alignment  with 
the  coordinate  system. 

(4)  Finally,  difficulties  in  the  stability  of  the 
results  may  arise,  especially  in  supersonic 
regions  in  the  finite-difference  method. 
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2.  BRIEF  DESCRIPTION  OF  THE  FLOW  DEVELOPMENT 
ABOUT  SWEPT  WING-BODY  COMBINATION 

Figure  1 shows  the  flow  development  about  a simple 
swept  wing-body  combination  with  increasing  free-stream 
Mach  number. 

At  moderate  subsonic  speeds,  if  the  lift  is  not  large, 
an  important  effect  of  the  sweep  is  to  move  the  loading 
forward  at  the  tip.  As  a result,  the  first  appearance  of 
supersonic  flow  occurs  in  the  tip  region.  Since  the 
isobars  near  the  wing  tip  lose  much  of  their  sweep,  the 
supersonic  flow  generally  gives  rise  to  a shock  wave-- 
the  initial  tip  shock  which  is  comparatively  weak  and  lies 
almost  normal  to  the  free  stream  [Figure  1(a)].  This  tip 
shock,  though  limited  in  the  spanwise  extent,  extends  to 
considerable  distances  above  and  below  the  wing. 

As  the  free-stream  Mro  is  increased,  the  initial  tip 
shock  moves  rearwards  over  the  wing  surface,  but  its 
influence  on  the  wing  is  limited  by  the  appearance  of  a 
second  shock..  This  rear  shock  rapidly  develops  to  affect 
a large  part  of  the  wing  span,  particularly  the  flow  ahead 
of  the  initial  tip  shock  [Figure  1(b)].  The  rear  shock  may 
be  regarded  as  associated  primarily  with  conditions  at  the 
wing  root.  It  consists  of  a compression  system  propagating 
outward  from  the  root  which  coalesces  on  the  outer  part  of 
the  wing  to  form  a shock  wave  [Figure  1(b)]. 

With  increasing  M^,  the  rear  shock  moves  aft  more 
rapidly  than  the  initial  tip  shock  which  is  overtaken  and 
disappears.  At  a sufficiently  high  Mw  the  rear  shock 
reaches  the  trailing  edge.  The  high  local  flow  velocities 
close  to  the  leading  edge  over  the  outer  part  of  the  wing 
lead  to  flow  separation  at  the  leading  edge.  This  starts 
near  the  tip  and  spreads  inboard  with  increasing  a.  For 
leading-edge  sweeps  greater  than  about  30  degrees,  the 
separated  flow  rolls  up  to  form  a partspan  vortex  lying 
obliquely  across  the  wing.  Above  a certain  stream  the 
flow  changes  in  type;  the  leading-edge  separation  is^ 
suppressed,  the  flow  passes  smoothly  around  the  leading 
edge  through  a forward  shock  wave  which  appears  to 
originate  close  to  the  leading  edge  at  some  spanwise 
position  [Figure  1(c)].  The  forward  shock  is  in  fact  the 
boundary  of  disturbances  from  the  inboard  part  of  the  wing 
leading  edge  which,  because  of  the  higher  local  supersonic 
Mach  number,  propagate  over  the  outer  wing  in  a direction 
more  highly  swept  than  the  leading  edge.  The  Mach  number 
component  normal  to  the  leading  edge  (M  cos  A)  for  which 
the  flow  attachment  occurs  varies  consi2erably , depending 
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on  A,  and  more  particularly  on  the  leading  edge  profile. 

In  general  M^  cos  A has  a value  between  0.55  and  0.90. 

As  the  free-stream  is  increased  beyond  that 
necessary  for  flow  attachment,  the  forward  shock  moves 
inboard  and  rearward  and  at  some  stage  intersects  the  rear 
shock.  Outboard  of  this  intersection  a strong  shock  forms 
(the  outboard  shock)  which  frequentlv  induces  a severe  flow 
separation  even  though  the  flow  is  attached  behind  both 
the  forward  and  rear  shock  [Figure  1(d)].  When  separation 
is  present  behind  either  of  these  shocks,  a vortex  tends 
to  form  and  the  associated  outflow  modifies  the  flow 
structure  outboard  and  becomes  an  important  factor  in  the 
ultimate  breakdown  in  the  attached  type  of  leading-edge 
flow  over  the  wing.  Separation  then  takes  place  along 
almost  the  entire  leading  edge,  and  the  shock  waves  no 
longer  have  a direct  influence  upon  the  surface  pressures 
which  are  dominated  by  a large  port-span  vortex  as  for 

lower  subsonic  M . 

00 

Many  of  the  features  described  above  are  shown  in  the 
surface-film  pattern  of  Figure  2,  obtained  at  M =0.95 
or  with  a = 10°. 

Complex  flow  patterns  similar  to  Figure  2 are  largely 
due  to  the  three-dimensional  nature  of  the  wing  flow  and 
hence  to  the  dominance  of  the  root  and  tip  influence  at 
transonic  speeds.  The  root  affects  the  flow  strongly 
behind  the  forward  shock  which  therefore  indicates  the 
limit  of  the  root  influence  and  the  flow  in  this  region 
is  partly  conical  in  character.  The  tip  influence  at 
transonic  speeds  is  delineated  by  small  disturbance  from 
near  the  tip  leading  edge  [the  tip  shock:  Figure  1(c)]. 
Ahead  of  the  forward,  outboard  and  tip  shocks  the  flow  is 
almost  two-dimensional  in  character.  The  extent  of  this 
zone  depends  on  the  shock  positions  which,  in  turn,  are 
dependent  on  wing  planform,  a,  and  M . By  making  the 
aspect  ratio  large,  or  by  deliberately  attempting  to  reduce 
the  root  and  tip  influence,  flow  can  be  obtained  over  a 
large  portion  of  the  wing.  This  flow  closely  resembles 
that  postulated  in  the  simple  sweepback  theory  which  is  so 
highly  desirable. 

The  relationship  between  swept-wing  flow  and  that  on 
the  equivalent  two-dimensional  section  at  a component  Mach 
number  normal  to  the  leading  edge,  M cos  A (if  not  too 
large)  is  shown  in  Figure  3.  At  higher  values  of  M cos  A 
it  is  difficult  to  minimize  the  effects  of  the  finite  wing 
aspect-ratio.  Recent  experiments  however  have  given 
satisfactory  results:  surface  pressures,  separation 


Equivalent  Local  Mach  Number  Component 
Normal  to  Leading  Edge 


M =0-50  M « 0-611  M -0  918 
a - 6°  a *4-91*  a *3  27* 


Figure  3.  Validity  of  Simple  Sweepback 
Concept  for  M cos  A = 0.5. 


boundaries  and  shock  positions  correlating  well  on  the 
basis  of  the  simple  sweep  theory. 


3.  DISCUSSION  OF  BASIC  EQUATIONS 


A.  Various  Forms  of  the  Transonic  Small  Perturbation 
Differential  Equations 

First  consider  the  exact  nonlinear  equation  for 
inviscid,  isentropic,  two-dimensional  potential  flow 
written  in  Cartesian  coordinates: 

(a2-4>2)4>  - 2$  <j>  $ + (a2-<J>2)$  = 0 

xx  x y xy  v yJ  yy 

where  a is  the  local  speed  of  sound,  and  $ , the  total 
velocity  potential.  In  terms  of  U = $x  and  V = $ , the 
total  velocity  components  of  Equation  (1)  become  y 

(a2-U2)Ux  = 2UV  Uy  - (a2  - V2)Vy  , 

V = U (irrotationality  condition) 
x y 

The  expressions  of  these  equations  in  terms  of  the 
perturbation  potential  $ and  the  perturbation  velocity 
components  u and  v,  defined  by 


U 1 + u ’ 


and  0 = x + <|>  (4) 

can  easily  be  written  down  in  full  (see  Reference  1, 
page  204).  We  will,  however,  consider  two  of  these 
perturbation  forms  for  our  guidance  in  selecting  the  best 
form  for  three-dimensional  flows. 

One  form  is  referred  to  as  the  small  perturbation 
equation  derived  by  neglecting  all  terms  containing 
products  of  perturbation  velocity  components.  We  will 
refer  to  this  as  NASA  Ames  ( transonic)  equation 
(Reference  2)  (NAE) 
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11-m2-(7+1)M^x]^xx  - 2M^y<t>xy  + [l-(Y-l)M^x]4,yy  = 0 (5) 


or  in  terms  of  u and  v: 

[l-M2-(Y+l)M2u]ux  = 2M 2 v uy  - [l-(Y-l)M2u]vy 

The  second  form  is  the  usual  classical  transonic 
equation  (CTE)  derived  from  Equations  (5)  and  (6)  by 
considering  the  special  behavior  of  the  fluid  flow  when 
the  local  velocity  is  close  to  the  speed  of  sound: 

[1-m2-(y+1)M^x]4>xx  + <f,yy  = 0 

or  [1-m2-(y+1)m2u]Ux  = - Vy 


v = u 
x y 


(6) 


(7) 


(8) 


We  will  now  analyze  these  equations  with  regard  to 
how  they  model  characteristic  lines,  sub-  and  supercritical 
flows,  and  shock  waves. 

When  Equations  (1)  through  (8)  are  used  for  supersonic 
flows,  they  can  be  analyzed  to  find  characteristic  lines; 
that  is,  lines  across  which  velocity  gradients  can  be 
discontinuous.  One  simple  way  to  derive  the  equation  for 
a characteristic  line  is  to  rotate  the  coordinates  through 
some  angle  0C  and  ask  if 
there  is  a value  of  0C  for 
which  an  advance  along  £ 
is  not  possible.  For 
Equation  (2) , it  is  easy 
to  show  that  a character- 
istic line  occurs  in  its 
solution  when 


/ 


/ 

/ 


9r 


[ —j  - cos  9 ] + [2v  cos0  sinO  -u(y-l  + 2cos  0 1 
^ c c c 

00 

2 2 

- [ (y-1) (U  ) + (u  cos0c-v  sin0c)2]  = 0 


(9) 


\ 


F 


where  the  brackets  enclose  terms  having  zero-,  first-  and 
second-order  powers  of  u and  v. 

Similarly,  it  can  be  shown  that  the  characteristic 
lines  of  Equations  (6)  and  (8)  are 

[- \ - cos20  ] + [2v  cos0  sin0  -u(y-l  + 2 cos20  )]  = 0 (10) 


and 

[\  - cos2d  ] - u (y+1)  cos 2 O = 0 (11) 

ir  c c 


respectively . 

Notice  that  Equations  (9)  and  (10)  are  identical  for 
both  the  zero  and  first  powers  of  u and  v.  But  Equation 
(11),  derived  from  the  classical  transonic  equation,  is 
not  in  agreement  with  Equation  (9)  even  through  the  first 
order  in  perturbation  velocities.  Since,  to  the  first- 
order,  oblique  shock  waves  bisect  characteristic  lines  of 
the  same  family,  this  result  can  influence  our  choice  of 
equations  when  we  seek  to  compute  flow  fields  with  oblique 
shocks . 


The  Critical  Velocity 

If  we  express  a second-order  partial  differential 
equation  in  the  form 

A $ +B0  +C<i>  + f(4>  , <J>  ,<t)  = 0 (12) 

xx  xy  yy  v x’  y*  ' v ' 

it  is  said  to  be  elliptic  or  hyperbolic,  depending  on 
whether  the  discriminant  (b2-4Ac)  is  <0  or  >0,  and  it 
changes  type  as  it  goes  through  zero.  This  occurs  in 
Equations  (1)  and  (2)  when  a2(Q2_a2)  goes  through  zero, 
where  |Q|  is  the  magnitude  of  the  local  velocity.  This 
gives  the  classical  result  that  Equations  (1)  and  (2) 
change  type  which  the  local  flow  speed  passes  through  the 
speed  of  sound  and  is  hyperbolic  on  the  supersonic  side. 
This  speed  is  referred  to  as  the  critical  speed  and  the 
value  of  the  pressure  coefficient 
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at  which  it  occurs,  as  the  critical  pressure  coefficient, 
C*. 

P 

Figure  4 shows  a plot  of  C*  versus  M„.  This 
relationship  is  valid  under  thepassumption  that  the  flow 
is  isentropic  along  a streamline  between  the  reference 
and  any  given  point.  Also  shown  in  Figure  4 is  the 
first-order  small  perturbation  approximation  to  = 2u*. 
This  approximation  follows  at  once  by  neglecting  all  terms 
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with  products  of  u and  v in  C -equation  and  can  be 
expressed  as  ^ 


- * 1 - €>2 


The  condition  under  which  Equation 

is  : 

m£(y2-1)4  - -'I2  4>x  Iy+1+(1-mJ)  (y-1)  ] + 


(14) 

(5)  changes  type 
(1-M2)  - M4  d)2  = 0 

N oo'  oo  T y 

(15) 


o 

There  are  two  roots  for  $*, 
physical  interest  leads  to 


and  the  one  of 


- 2u* 


M 


-5 


Y-l/i  \<2 \ 

7+r(1_*v 


] 


(16) 


2 

Equation  (16)  reduces  to  Equation  (14)  if  A and  higher 


order  terms  in  p are  neglected. 


It  is  easy  to  show  that  the  expression  for  -2u* 
derived  from  Equation  (7)  is  identical  to  Equation  (14). 
We  therefore  conclude  that,  to  the  lowest  order  in 
perturbation  velocities,  Equations  (1),  (2),  (5),  (6), 

(7)  and  (8)  all  change  type  according  to  the  same 
relationship  among  u,  y and  as  given  by  Equation  (14). 


Simple  Sweep  Theory 

The  simple  sweep  theory  is  generally  derived  by 
considering  a sheared  wing  of  infinite  span  and  constant 
section  (Reference  3) . The  principal  assumption  is  that 
the  velocity  component  parallel  to  an  edge  is  constant  or 
that  the  perturbation  velocity  parallel  to  an  edge  is  zero. 
Under  these  conditions  a shock,  if  it  occurs,  would  also 
have  to  be  parallel  to  an  edge,  and  hence  its  strength 
would  depend  only  on  the  component  of  velocity  normal  to 
an  edge.  We  will  use  the  symbol  C*(0)  to  designate  the 
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pressure  coefficient  in  an  isentropic  flow  at  which  the 
component  of  the  total  velocity  normal  to  some  plane,  fixed 
by  6 = constant,  is  sonic.  This  plane  will  be  assumed  to 
be  vertical,  i.e.,  perpendicular  to  the  xy-plane.  Note 
that  Cp(0)  = Cp , which  has  the  same  meaning  as  given  in  the 
previous  section.  " 

If  then  q(0)  is  the 
component  of  velocity 
normal  to  an  isobar  plane, 
classical  sweep  theory 
leads  at  once  to  the 
critical  value  of  q(0) : 


q*2(0)  - U2  sin2e  = a2 


(17) 


From  Equation  (14)  it  can  be  shown  that,  for  isentropic 
flow 


1 . , 


2 cos  0 , 1 




W 


l) 


M cos  0 


(18) 


Equation  (18)  can  be  used  to  derive  the  exact  value  for 

C*(0): 


yM‘ 


c*(0)  = [(4r  (1  + ^ M2  cos20)}y 


y^T 


1] 


(19) 


Applying  the  small  perturbation  approximation  to 
Equation  (19)  yields 


- 2u*(9)  * 1 - [ 


q*(0) 


2 cos  0 / 1 

— 


7+T 


Y~) 

cosz0 


(20) 


Figure  5 compares  Equation  (19)  with  Equation  (20) , 
which  is  the  simple  sweep-theory  counterpart  of 
Equation  (14) . 
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Figure  5.  Cp(9)  versus  0 as  Given 

by  Equations  (19)  and  (20) . 


Shocks 

So  far  the  discussion  has  been  restricted  to  the 
properties  of  genuine  solutions  of  the  differential 
equation,  i.e.,  solutions  that  result  in  continuous 
variations  in  the  velocity  components.  We  must  next 
inspect  the  properties  of  possible  weak  solutions  of 
these  equations  (see,  for  example,  Reference  4).  A weak 
solution  can  be  composed  of  two  genuine  solutions 


separated  by  a surface  across  which  the  dependent  variables 
can  be  discontinuous.  We  refer  to  this  surface  as  a shock, 
regardless  of  the  form  of  the  hyperbolic  equation  being 
investigated.  For  the  Eulerian  equations,  the  jump  condi- 
tions for  the  dependent  variables  across  the  shock  are 
referred  to  as  the  Rankine-Hugoniot  conditions.  In 
general,  the  conditions  that  must  be  satisfied  across  a 
shock  can  be  found  by  writing  the  partial  differential 
equations  in  conservation  form  and  connecting  the 
differences  in  the  conserved  variables  with  the  direction 
cosines  of  the  surface  across  which  they  are  discontinuous 
(Reference  4).  The  jump  conditions  for  the  isentropic 
shock  embedded  in  Equations  (1)  and  (2)  lead  to  an  implicit 
relationship  between  the  velocity  components  on  the  two 
sides  of  the  discontinuity,  which  can  be  written  as 
(Reference  5) : 

Of  Ul  - G2  U2  = " G2v2^tan0  (21) 

where 

1 

y-1  2 2 

Gi  = 11  - V:(ui+vi)]  : i = 1.2 

and  the  subscripts  1 and  2 refer  to  conditions  on  the 
upstream  and  downstream  sides  of  the  shock  at  any  given 
point.  The  local  slope  of  the  shock  is  -tan0. 

The  conservative  form  of  Equation  (8)  can  be  written 
as 


mV  - \ M2  (y+1)u2]  + ^ tv-M2  (y-l)uv]  = 0 

(22) 


condition,  we  obtain 
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^ [(1-M2)u-^1 


9v  3u  _ n 
9x  ” 3y  u 


and  for  its  jump 


[1  - ^ M2  (y+1)  (u^+u2)  ] tan20  - 2 M2 (v2+u2tan0) 


+ 1 ~ - 7 m2(y+1)(Ui+u2)  = 0 (23) 


The  jump  condition  for  the  classical  transonic 
equation,  [(7)  and  (3)]  is  important  for  our  purpose. 
The  conservative  form  of  Equations  (7)  and  (8)  is 


8 

3x 


[(l-Ml)u  - 


l£(Y+l)u2] 


3v 

9y 


= 0 


(24) 


3 v _ 3u 
3x  3y 


and  the  jump  conditions  yield  the  equations: 


(l-M2)Au  - j M2(y+1)Au2  = Av  tan0 


Av  = - Au  tan0 


where 


Au  = u-^  - u2  ; Av  = v-^  - v2 

Eliminating  Av,  by  noting  that 

2 2 2 — 

A(u  ) = u^  - u2  = (u^-u2) (u^+u2)  = (u^+u2)Au  = 2u  Au 

we  have 

tan20  = M2  - 1 + M2(y+1)u  (25) 


or 


(-2u1)+(-2u2) 

2 


00 


1 

T~ 

cos  9 


1) 


(26) 


\ 
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This  is  the  condition  that  must  be  met  across  the  shock 
according  to  the  classical  transonic  theory  expressed  by 
Equations  (7)  and  (8). 


Shocks  and  Supercritical  Flow 

A flow  becomes  supercritical  when  the  differential 
equation  changes  type  from  elliptic  to  hyperbolic.  To 
the  lowest  order  in  the  perturbation  velocities,  it  has 
been  shown  that  this  occurs  under  the  same  condition  for 
Equations  (1)  through  (8),  namely,  where 

- 2u  = - 2u*  - - (*T  - 1)  (27) 

Y+1  Mz 


When  all  simulated  shocks  are  nearly  perpendicular 
to  the  free-stream  direction,  that  is  to  say  the 
approximation  9=0  and  = 0 is  good  enough,  Equations 
(22)  and  (24)  lead  to  thezsame  result.  This  result  applies 
to  the  standard  two-dimensional  flow  simulations  and,  in 
the  form  of  Equation  (26) , it  gives  a condition  that  must 
be  met  by  the  average  of  the  perturbation  velocities 
across  the  discontinuity.  Notice  that  as  u,  -*■  u? , they 
both  approach  u*.  This  means  that  the  intensity^of  a 
shock  perpendicular  to  the  free  stream  is  centered  about 
the  critical  perturbation  velocity,  which,  in  turn,  is 
based  on  the  condition  that  the  equation  changes  type. 

For  such  cases,  Equations  (5)  through  (8)  give  essentially 
the  same  results,  and  Equations  (7)  and  (8)  are  to  be 
preferred  because  of  simplicity. 

If  the  boundary  conditions  are  such  that  a vertical 
shock  is  oblique  to  the  free  stream,  the  shock  position 
and  intensity  could  be  significantly  different  depending 
on  whether  Equations  (5)  and  (6)  or  Equations  (7)  and  (8) 
were  used  for  the  simulation.  For  example,  consider  the 
special  case  where  v2  = - 2u2  tanG , which  exists  when  the 
component  of  the  perturbation  velocity  parallel  to  the 
shock  is  zero,  which  is  the  condition  for  simple  sweep 
theory.  Under  this  condition,  the  relation  [Equation  (23)] 
reduces  to 


(-2ux)+(-2u2) 

2 


2 cos^G 

~Tfl — 


1 

2~ 

cos  9 


1) 


(28) 


\ 
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On  the  other  hand,  the  jump  condition  for  Equations 
(7)  and  (8) , given  by  Equation  (26)  , does  not  depend  on 
any  special  relation  between  U2  and  V2  and  is  therefore 
valid  for  all  values  of  0.  This  simple  treatment  brings 
out  clearly  the  difference  between  Equations  (5)  and  (6) 
and  Equations  (7)  and  (8)  in  simulating  oblique  shocks 


Figure  6.  Exact  and  Approximate  Values 
of  Cp(0)  for  Mro  = 0.85. 
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from  a comparison  of  Equations  (28)  and  (26) . This 
difference  is  shown  in  Figure  6 for  M = 0.85.  Exact 
value  of  Cp (9)  from  simple  sweep  theory,  Equation  (19), 
is  also  plotted  in  Figure  6.  It  is  clearly  seen  from 
Figure  6 that  Equations  (7)  and  (8)  are  very  poor  models 
for  flows  with  shocks  that  are  more  than  30  degrees 
oblique  to  the  free  stream.  Equations  (5)  and  (6),  on 
the  other  hand,  are  acceptable  for  a wide  range  of  0, 
limited  only  by  the  accuracy  of  the  small  perturbation 
approximation  itself. 

Hence,  for  three-dimensional  transonic  flows  in  which 
the  presence  of  oblique  shocks  at  moderate  to  large  sweep 
angles  is  anticipated,  the  proper  form  of  the  transonic 
small  disturbance  equation  is 

| 1 1-M l-  (y+l)M^x]  4>xx  - 2M^y<Pxy  + 1 1-  (Y-1)M20)X]  4>yy  + 4»zz  = 0 

(29) 


SECTION  II 


THREE-DIMENSIONAL  TRANSONIC  SMALL-PERTURBATION  EQUATIONS 
AND  THE  ORDER  OF  MAGNITUDE  ANALYSIS 


Consider  first  the  exact  equation  for  the  velocity 
potential  $ in  rectangular  coordinates  x' , the  local  stream 
direction  and  y'  and  z'  normal  to  the  free  stream: 


(a  "9  ) ^x'x'+  a % ' y ' a ®z'z’  — ^ 


(30) 


where  a and  q are  the  local  values  of  the  speed  of  sound 
and  stream  speed,  respectively,  related  by 


X^i  q2 


(31) 


where  Mra  is  the  free-stream  Mach  number  and  the  froe-stream 
velocity  has  been  taken  to  be  of  unit  magnitude.  The 
underlined  term  in  Equation  (30)  determines  when  the 
equation  changes  type. 

A transformation  from  the  (x1 ,y' , z ') -system  to  the 
system  (x,y,z)  is  now  made  by  rotating  the  coordinate  axes. 
Here  x is  the  free-stream  direction  and  y and  z normal  to 
the  free  stream: 


~2(a2-q2)  [u2$xx  +v2*yy  +w2*zZ 


+ 2uv»xy  + 2uwlxz 


+ 2vw$yz] 


+ ^ [ (q2‘U2)<txx  + (q2-v2U  + (q2-w2)$zz 

q 

- 2uv$xy  - 2uw<I>xz  - 2vw$yz  ] =0  (32) 


where^u  = <I>X,  v = $y  and  w = $z  are  the  velocity  components 
and  q = u2  + v2  + w2 . In  Equation  (32)  the  underlined 
terms  result  from  the  underlined  term  in  Equation  (30) . 

This  equation  is  no  more  than  a rearrangement  of  the 
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standard  compressible  flow  equation  for  <t> , in  which  each 
of  the  second  derivatives  has  been  split  into  two  parts. 


Transonic  Small-Perturbation  Approximation 

The  transonic  small-perturbation  equation  may  be 
derived  from  the  exact  Equation  (32)  in  a number  of  ways. 
We  give  here  a derivation  in  which  the  scaling  of  the 
variables  is  avoided  but  the  physical  meaning  of  the  basic 
assumptions  is  stressed. 

The  total  velocity  potential  $ is  first  replaced  by 
a perturbation  potential  <|>  defined  by 


$>  = x + 


(33) 


where  = 1.  The  velocity  components  are  now 

u = 1 + <J>X  , v = 0 , w = 4>z  . (34) 

We  assume  first  that  <J>X,  <f>y  and  <p  are  all  small  and 
if  the  z-axis  is  taken  in  the  verticalairection  the  slope 
of  the  wing-body  surface  relative  to  the  horizontal  plane 

z = constant  is  approximately  given  by  £-5-  = — . We  shall 

W oX  U 

assume  that  — is  small: 


w 

u 


6 


6 <<  1 


(35) 


This  small-perturbation  approximation  cannot  obviously 
be  valid  near  a blunt  leading  edge  but,  as  in  two- 
dimensional  sections,  we  accept  this  local  inconsistency 
and  proceed.  It  is  convenient  to  choose  the  streamwise 
chord  length  to  be  of  unit  length.  So  the  condition 
6 <<  1 implies  that  <p  <<  1 and  we  write 

4>  ~ e , e « 1 . (36) 

Now  we  introduce  two  lengths  b and  t in  the  y and  z 
directions,  respectively,  such  that 
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The  length  b is  related  to  the  planform  of  the  wing.  Thus 
for  any  wing  of  aspect  ratio  ~1  or  of  any  aspect  ratio  but 
appreciable  sweep,  b ~ 1.  On  the  other  hand,  for  a wing 
of  high  aspect  ratio  and  small  sweep,  b >>  1.  This  is 
consistent  with  the  exact  relation  <py  = - <px  tanA  for  an 
infinite  swept  wing  at  an  angle  A with  constant  chord 

length.  The  length  t is  related  to  the  wing  slope  jj  ~ 6 , 

since  4>x  = ^ approximately,  so  that 


£ 

t 


~ 6 


(38) 


The  next  step  is  to  expand  the  terms  in  Equation  (32) 
in  the  new  dependent  variable  <j) , by  making  use  of  the 
conditions : 


4>  ~ e <<  1 


1 


1 


(39) 


Thus,  on  omitting  quantities  of  third  and  higher  order,  we 
have  Equations  (31) , (32) , (34)  and  (39) : 


[1-M2- (y+l)M24>  ] 4>  + 2 (1-M2)  [ d>  <t>  +<b  d>  ] 

1 oo  v 1 7 °°rx  xx  v °°  LYyYxy  Yz  xz 

+ [ 1-  (y-1)M2(J>  ] (<p  +<p  ) - 2<p  <fi  - 2d)  4>  =0 

1 VI  ' ooYx  v yy  zz7  y xy  Yz  xz 

(40) 


To  the  first-order,  the  pressure  coefficient  is  given  by 

Cp  = - 2*x  (41) 

The  small-perturbation  Equation  (40)  is  now  further 
simplified  by  making  a transonic  flow  approximation.  We 
suppose  that  1 - m£  «,  $ that  is 


21 


1 


(42) 


e 


Thus  the  second  of  the  underlined  terms  in  Equation  (40) 

o 

can  be  omitted.  Furthermore  the  term  (Y“l)^00<))x(4>yy+4IZ2) ' 

can  also  be  omitted.  The  result  is  the  transonic  small- 
disturbance  equation  for  arbitrary  wing-body  combination: 


[ l-M2- (v+1)M2»xhxx  + »vy  + - 2 Vxy  - 2«z»xz  - 0 (43) 


This  equation  differs  from  the  commonly  adopted  equation 

in  having  additional  terms  2<j>  . <j>  and  2<p 

y xy  z xz 


Note  that  for  b ~ 1,  that  is,  for  wings  of  aspect 
ratio  1,  or  for  wings  of  any  aspect-ratio  but  appreciable 


sweep,  (<j>yy  + <(>zz)  << 


Yet  <t> 


that  <p 


yy 


~ e 


so  it  follows 


zz 


£ and  the  length  t ~ 1.  However,  for  b 


that  is,  for  two-dimensional  flows 

A / 


while  d>  <b 

zTxz 


(~~t)  <<: 


0, 


<b  -*■  0 and  6 d> 
yyy  ^y^xy 

. „„(~-V)  , and  the  equation  reduces  to 
z.  z z z 

the  well-known  two-dimensional  nonlinear  transonic  form: 


[1-M2-(y+1)M *<Px]<Pxx  + 4>xx  * 0 (44) 

which  implies  that  $ e2  so  that  t ~ e”1^2,  that  is 

t >>  1.  From  the  relations  (36)  and  (38)  it  follows  that 
4>  — 1 5 and  is  larger  for  a two-dimensional  airfoil  than 
for  a finite  wing  of  the  same  thickness. 

The  magnitudes  of  various  quantities  for  three-  and 
two-dimensional  wings  follow. 
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MAGNITUDES  OF  VARIOUS  QUANTITIES  FOR 
THREE-  AND  TWO- DIVISIONAL  WINGS 


Quantity 

3-Dim.  Wing 

2-Dim.  Wing 

b:  transverse  y-direction 

length  scale 

1 

00 

t:  length  scale  in  z-direction 

1 

6-1/3 

<f> : perturbation  potential 

6 

62  / 3 

1-M2 

OO 

6 

62/3 

4>  , ~ pressure  coefficient 

6 

6 2 / 3 

d>  , streamwise  pressure 

gradient 

<t>  , vertical  pressure 

gradient 

6 

6 2 / 3 

6 

6 

Mxx 

62 

64/3 

( d>  +d>  ) 

yyy  Yxz; 

terms  in  governing 

r 2 

r 4/3 

differential 

0 

0 

^Z^XZ  j 

equations 

62 

82 

Comments 


For  the  two-dimensional  airfoil,  the  magnitudes  that 
have  emerged  for  t,  <})  and  l-M^  show  that  if  new  scaled 
quantities  of  unit  order  of  magnitude  are  to  be  defined 
then  the  appropriate  scaling  is  obtained  by  writing 


K 


a'2/3 

1-M2 

00 

7m 


(45) 
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Equation  (44)  then  reduces  to 

[K  - (y+1HxHxx  + 4-z£  = 0 

For  a three-dimensional  wing  the  corresponding 
scaled  quantities  are  different.  They  are 

z ~ z 

<F  - 6_1<|) 


(46) 


(47) 


It  is  very  important  to  recognize  these  differences. 
Note  that  from  inspection  of  Equation  (43)  it  does  not 
seem  possible  to  define  similarity  parameters  like  K such 
that  the  solution  depends  only  on  K (as  with  two-dimensional 
airfoils)  and  not  on  the  thickness  6 as  well.  In  other 
words,  there  does  not  appear  to  be  a simple  transonic 
similarity  rule  for  finite  wings  as  there  is  for  airfoils. 

The  differences  in  the  vertical  length  scale  t and 
the  perturbation  potential  <(>  imply  differences  in  the 
physical  flow  fields.  The  pressures  and  streamwise 
pressure  gradients  on  the  wing  surface  are  smaller  by  a 
factor  of  ~ §1/3  than  those  on  the  corresponding  airfoil 
surface.  For  a finite  wing  the  normal  and  streamwise 
pressure  gradients  are  of  the  same  order  of  magnitude,  but 
for  an  airfoil  section  the  normal  gradients  are  smaller  than 
the  streamwise  gradients  by  a factor  of  ~ §1/3  . 
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V 


SECTION  III 


SUITABILITY  OF  THE  VARIOUS  TRANSONIC  EQUATIONS 


We  have  already  derived  the  NASA  Ames  equation 
given  by  Lomax,  Bailey  and  Ballhaus  (Reference  2) : 

[(1-I-6-(y+1)M^x1*xx+  (I-(y-1)M^x  »yy  + *zz-2^ty»xy  - 0 

(43) 

This  equation  was  suggested  for  use  in  the  case  of  flow 
over  finite  wing-body  combinations  where  oblique  shocks  at 
moderate  to  large  sweep  angles  are  anticipated.  If  6 is 
the  thickness-to-chord  ratio,  the  assumptions  leading  to 
this  equation  are : 


4>  = 0 ( 6 2 7 3 ) , 
z-scaling  according  to  z = 
and  (1  - M^)  = 0(1)  . 

For  flow  with  shocks  we  have 

M2cos29-l  = M2cos20-1  + 2M2(1  +^M2)<J>  cos20  + . . . > 0 

CO  co  / oo  ' T v 


(49) 


(50) 


where  M is  the  local  Mach  number  in  front  of  the  shock  and 
9 the  shock  sweep  angle.  Now  we  have  4>x  = 0(62/3)  > o and 
consequently  for  Mot  < 1 Equation  (50)  requires 

M2  cos20  - 1 = 0(6P)  , p > 2/3 


or 


M 


2 - 1 = tan20  + 0(6p)cos29 


p > 2/3 


(51) 
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whence 


Mjj  - 1 - M2  cos2e  - 1 = 0(62/3)  . 


The  ignored  entropy  rise  across  a shock  is  then  of  order 
(^n" ^ ^ ~ = 62  • which  is  consistent  with  the  fact 

that  Equation  (48)  is  a second  order  accurate  approximation. 


Now  since  we  expect  normal  shocks  to  occur  in  any 
case , we  may  have  0=0  and  thus 

M2  - 1 = O(fiP)  , p * 2/3  , 


O 

which  contradicts  the  assumption  1 - = 0(1) . This  means 

that  Equation  (48)  can  only  be  used  for  subcritical  flow; 
that  is  to  low  transonic  speed  regime.  Yet,  it  is  a 
nonlinear  equation  and  as  such  capable  of  producing  shocks. 

The  classical  Guderley-Karman  equation  is 


[1-M2- (y+l)M2<f>  ]d>  + d>  + d>  = 0 

L co  ' ' «rxjyxx  Tyy  vzz 


(52) 


which  is  usually  given  for  a slender  body  as  well  as  for 
high  aspect-ratio  wings  of  very  small  sweep  and  dihedral. 
The  necessary  assumptions  are 

(!)  = 0(62/3)  ; 


; 


y-  and  z-scaling  is  according  to 


and 

1 - M2  = 0(62/3) 

whence  we  attribute  it  to  the  medium  transonic  speed  regime 
since  M*,  is  closer  to  one  than  in  the  foregoing  case.  In 
the  derivation  of  this  equation  no  contradictions  are 
encountered  for  the  case  of  flow  with  shocks  as  long  as 
they  are  slightly  swept. 
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Newman  and  Klunker  (Reference  6)  have  added  one 
third  order  term  to  the  Guderley-Karman  equation  for  a 
better  approximation  of  the  critical  speed  where  the 
equation  changes  type  from  elliptic  to  hyperbolic: 

1(1-m£)-(y+1)!&x-  \ (y+1)M^1  txx+tyy  + 4>zz  = 0 (54) 

The  necessary  assumptions  are  the  same  as  given  in 
Equation  (53) . 

In  Reference  7 (see  also  Reference  8) , Hall  and 
Firmin  have  derived  and  used  the  equation  (see  Section  4) : 

[<l-ME)-<Y+l)Hi*x]»xx  + *yy  + *2z  - 2,ytxy  - 2*2tX2-0  (55) 

for  a low  aspect-ratio  wing  fuselage  combination  or  a 
wing  of  any  aspect-ratio  but  appreciable  sweep.  The 
necessary  assumptions  are 

<)>  = 0(6)  , 1 - = 0(6)  (56) 

whence  we  attribute  it  to  the  high  transonic  speed  regime 
since  Moo  is  still  closer  to  one  than  in  the  foregoing 
case.  Also  in  the  derivation  of  this  equation  no  contra- 
dictions are  encountered  for  the  case  of  flow  with  shocks. 


The  Conservation  Forms  of  Small  Perturbation  Transonic 


The  mass  conservation  law  for  the  compressible  fluid 
flow  may  be  written  exactly  as: 

[ p (l+4>x)  ] + ( p<(>  ] + [ P 4> z ] = 0 (57) 


Various  approximations  to  Equation  (57)  may  be 
derived  by  substituting  the  properly  truncated  series  for 
the  terms  inside  the  brackets: 
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°<l+*x>  - 1 + U-niHx  - ^3-(2-y)M^M^2 


7 " 7 M2<t>2  + 

2 00  y 2 °oyz 


(58) 


ry  " " ^x't’y  + 


P<t>  = <p  - M (p  <p  + . . 
z rz  00  xYz 


in  Equation  (57).  For  instance,  keeping  only  terms  up  to 
the  second  powers  of  <p  — derivatives  in  the  x-  and  y-directions 
the  conservation  form  of  the  transonic  equation  given  by 
van  der  Vooren  and  his  associates  at  NLR  in  Holland 

[a-t64.x-^3-(2-Y)M2}M^2-lM^2]  + Uy-M£+X*y]  + [+ZJ  -0 

x y z 

(59) 


The  conservation  form  of  the  NASA  Ames  Equation  (48) 
turns  out  to  be 


Ul~M *)<Px  - ^(Y+1)M^  + £(y-3)M^] 

y x 

+ [<i>  -(y-Dm^a  ] + [cj>  ] = o (60) 

3 y z 

The  Equation  (59)  can  be  shown  to  differ  from 
Equation  (60)  only  by  third  order  terms.  However  their  shock 
relations  differ  significantly  and  will  be  compared  below. 


Investigations  of  Shock  Relations 

Consider  three  typical  cases  of  vertical  shocks,  for 
which  <pz  = 0 . 

The  first  case  is  the  normal  shock  which  occurs,  for 
example,  at  the  wing  roots.  We  have  <)>  = <p  = 0 . The 

y 2 
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results  are  shown  in  Figure  7 for  0.8  1 1 1.  The  NASA 

Ames,  Guderley-Karman  and  RAE  equations  [Figure  7(a)]  show 
fair  agreement  with  full  potential  equation. 


a)  NASA  AMES,  GUDERLEY- VON 
KARMAN.  RAE  f~0  . 

(p-O.  00  ^M.pSD 


C)  Kl.UNKER-NF.WMAN  (*) 
(p.O,  00  5iM„  =51) 


(-20*  =5  p =5*20*.  M^.064) 


Figure  7.  Normal  Shock  Relations 
for  Various  Transonic  Equations. 
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Figure  8(a).  Infinite-Sheared  Shock 
Relations  for  Various  Transonic  Small 
Perturbation  Equations  (Mw  = 0.84,  A = 30°)  . 
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Figure  8(b).  Shock  Relations  for  Various 
Transonic  Small  Perturbation  Equations  in 
a Case  of  Forward  Shock  (Moo  = 0.84,  A = 30  ). 
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Figure  8(c).  Upper  Surface  Pressure  on  C141 
Panelmodel  for  Various  Transonic  Small 
Perturbation  Equations. 
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Figure  8(d).  Pressure  Distribution  on 
ONERA  M6  Wing  Model. 
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!For  values  of  $ 1.3  there  is  also  fair  agreement 
with  the  Rankine-Hugoniot  relations.  The  results  of  NLR 
equation  [Figure  7(b)]  are  practically  independent  of  M^, 
and  agree  well  with  the  potential  theory.  The  results  of 
tne  Klunker-Newman  equation  [Figure  7(c)]  are  nearly 
independent  of  Moo.  Though  the  agreement  with  full  potential 
theory  as  well  as  with  the  Rankine-Hugoniot  relations  is 
good  for  low  M^-values  the  overall  picture  is  not  as  good 
as  for  the  NLR  equation,  or  even  the  other  three  equations 
for  reasons  already  indicated. 

Secondly,  the  strong  oblique  shock  on  an  infinite 
sheared  wing,  for  the  various  equations,  was  compared. 

If  A is  the  sweep  angle  of  both  wing  and  shock  we  have 
<t>y  = " 4>x  tan  A.  The  results  of  the  specific  case 
Moo  = 0.84,  A = 30°  are  shown  in  Figure  8(a).  It  is  apparent 
that  neither  the  Ouderley-Karman  nor  the  Klunker-Newman 
equation  can  describe  this  kind  of  shock  adequately.  The 
results  obtained  with  NASA  Ames  and  RAE  equations  are 
reasonable,  the  RAE  equation  being  slightly  better.  The 
NLR  equation  gives  good  agreement  with  full  potential  theory 
and  for  M]_  £ 1.3  fair  agreement  with  the  Rankine-Hugoniot 
relations.  For  Ma>  = 1 , the  NASA  Ames  and  RAE  equations 
also  agree  very  well  with  full  potential  theory. 

The  third  case  of  interest  is  the  weak  oblique  forward 
shock.  Assuming  infinite-sheared  conditions  (<f>„  = - <px  tan  A) 
in  front  of  the  shock  and  flow  parallel  to  the  fuselage 
(4>y  = 0)  behind  the  shock.  The  results  for  the  case 
Mco  = 0.84,  A = 30°  are  shown  in  Figure  8(b).  The 
conclusions  are  the  same  as  for  the  second  case.  Only  the 
NLR  equation  gives  fair  agreement  with  both  the  full 
potential  equation  and  the  Rankine-Hugoniot  relations. 

In  Reference  9,  van  der  Vooren  and  his  associates  have 
given  some  results  of  calculations  of  transonic  flows  over 
isolated  semi-wings.  Nonconservative  as  well  as  fully - 
conservative  rotated  difference  schemes  have  been  used. 

They  came  to  the  conclusion  that  a nonconservative  scheme 
is  inadequate  for  shock  capturing  capability  [see  Figures 
8(c)  and  8 (d) ] . 
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SECTION  IV 

INTEGRAL  EQUATION  FOR  THREE-DIMENSIONAL  TRANSONIC 
FLOWS  PAST  WING -BODY  COMBINATION 

THE  NONLIFTING  CASE 


Theoretical  pressure  distributions  on  nonlifting 
circular  arc  airfoils  immersed  in  a high  subsonic  free 
stream  were  first  given  by  Oswatitsch  in  1950  by  approxi- 
mately solving  the  integral  equation,  which  was  proposed 
by  him.  The  analysis  is  carried  out  in  the  physical  rather 
than  in  the  hodograph  variables,  and  leads  to  a nonlinear 
integral  equation  in  which  the  unknown  velocity  appears 
outside  as  well  as  inside  the  integral.  Oswatitsch  found 
approximate  solutions  not  by  iteration,  but  by  introducing 
various  functions  containing  undetermined  parameters  into 
the  integral  equation  and  then  by  determining  these 
unknowns  by  satisfying  the  integral  equation  at  a small 
number  of  points  on  the  airfoil  chord.  The  application  of 
the  method  to  some  airfoil  section  did  show  certain 
definite  characteristics  of  transonic  flow  such  as  the 
appearance  of  shock  waves  and  their  rearward  movement 
across  the  chord  with  increasing  Mach  number.  However,  the 
method  fails  to  give  proper  results  at  high  subsonic  M«, 
greater  than  about  0.88  with  a 6 percent- thick  circular 
arc  airfoil. 

Gullstrand,  at  KTH,  Sweden,  tried  to  rectify  the 
situation  by  seeking  a solution  by  iteration.  Gullstrand 
uses  the  integral  equation  to  determine  only  the  solution 
for  the  forward  part  of  the  airfoil  and  then  uses  the  method 
of  characteristics  to  complete  the  solution  for  the  rear  of 
the  airfoil. 

Spreiter  and  Alksne  (Reference  10)  proposed  in  1955 
a solution  of  the  three-dimensional  integral  equation  given 
by  Oswatitsch  using  an  iteration  process  permitting  the 
integral  equation  to  be  satisfied  at  a much  larger  number 
of  points  than  in  the  original  method  of  Oswatitsch.  This 
method  gives  approximate  solutions  at  all  Mach  numbers  up 
to  unity,  and  appears  to  avoid  any  multiplicity  of  solutions 
that  were  unavoidable  in  Oswatitsch' s earlier  method  for 
supercritical  Mach  numbers. 
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Basic  Equations 

The  simplified  small  perturbation  equation  in  the 
transonic  flow  regime  is  the  usual  nonlinear  equation 

+ *yy  + *zz  - K *X*XX  <61> 


where 


k 


H 


2 y+1 
U 


(62) 


and  (p  , the  perturbation  velocity  potential  given  by 


<p  = - U^x  + $ 

As  a result 
of  minor  variations 
in  the  perturbation 
analysis,  other 
authors  have  used 
at  least  four 
different  relations 
for  k,  the  coeffi- 
cient of  the  non- 
linear term  in  the 
simplified  equation 
for  the  transonic 
flow.  A straight- 
forward development 
of  the  second-order 
theory  leads  to  the 
relation  given  in 
Equation  (62) . 

This  is  sometimes 
simplified  to 


Figure  9.  View  of  Wing 
and  the  Coordinate 
System. 


k = 


Y+1 


U 


CO 


(63) 


by  arguing  that  Moo  can  be  set  equal  to  unity  since  the 
right-hand  side  is  after  all  an  approximation  to  the 
actual  second-order  equation.  If  the  full  nonlinear 
equation  in  i>  , 
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(a2-u2)<}>  - 2uv<{>  + (a2-v2)<t 

x 'Yxx  Yxy  yy 

2 

the  total  velocity  potential  is  divided  by  a and  the 
quotient  l/az  in  each  term  is  expanded  in  a binomial 
series,  the  resulting  coefficient  k of  the  term  <f'x$xx 
turns  out  to  be 


k = ll2 

oo 


'2  + (y-1)m; 
u~~ 


2n 


(64) 


Still  another  expression  for  k is  used  by 
Oswatitsch  by  writing 


i-m: 


a*-lT 


+ . 


thereby  giving  the  value  of  k as 

,2 


k = 


1-M 

OU 

aW-V 


(65) 


where  a*  is  the  critical  sound  speed  given  by 


— = ■ /y~1  i 2 
uoo  Y Y+T  M2(y+1) 

A similar  situation  arises  in  the  derivation  of  the 
simplified  equation  for  the  shock  polar  given  below.  A 
significant  case  where  the  four  relations  for  k lead  to 
different  results  is  the  prediction  of  the  variation  of  the 
critical  pressure  coefficient  Cpcr  with  M^.  The  critical 

pressure  coefficient  Cpcr  is  the  value  of  Cp  at  a point 

where  local  M = 1.  This  condition  is  recognized  by  the 
vanishing  of  the  coefficient  of  4>xx,  thus 
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1 - M - k(4>  ) = 0 

oo  vyxy 
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or 
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(66) 


The  exact  relation  for  isentropic  flow  from 
Reference  1 is 


C 


P 


cr 


M 


1 


Similar  comparisons  can  be  made  for  local  Mach 
numbers  M other  than  unity  or  by  considering  the  velocity 
jump  across  a shock  wave  or  by  comparison  to  some  of  the 
existing  experimental  results  on  drag  of  a wedge.  Each 
of  these  comparisons  provides  striking  evidence  in  support 
of  the  value  of  k given  by  Equation  (62)  although  in  some 
cases  Oswatitsch's  value,  Equation  (65),  gave  good  results 
as  well. 

Equation  (61)  is  valid  only  in  regions  where  the 
necessary  derivations  exist  and  are  continuous.  Since 
these  conditions  do  not  hold  where  shock  waves  occur,  an 
additional  equation  is  needed  for  the  transition  through 
the  shock.  The  necessary  equation  is  provided  by  the 
classical  regulation  for  the  shock  polar: 


~2  ~2 
v^  + w^ 


= (u1-u2)‘ 


ulu2 


- a* 


u. 


uiu2  + a* 


(67) 


where  u,  v and  w refer  to  Cartesian  velocity  components 
with  u being  parallel  to  the  flow  direction  ahead  of  the 
shock,  the  subscripts  1 and  2 refer  to  conditions  ahead  of 
and  behind  the  shock,  and  a*  is  the  critical  sound  speed. 
Carrying  out  a small  perturbation  analysis  of  Equation  (67) 
analogous  to  that  performed  in  the  derivation  of  Equation 
(61) , it  can  be  shown  that  the  following  relation  results 
between  the  perturbation  velocity  components  on  the  two 
sides  of  the  shock  wave : 
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U-M*)  (uru2)2  + (v1-v2)2  + (wL-w2)2  = k 


ul+u2  , ,2 

~T~  (ul-u2> 


where  now  u,  v and  w are  the  perturbation  velocity 
components  parallel  to  x,  y and  z-axes.  This  equation 
corresponds  to  the  shock-polar  curve  for  shock  saves  of 
small  strength  inclined  at  any  angle  between  that  of 
normal  shock  waves  and  that  of  the  Mach  lines.  On  either 
side  of  the  shock  wave,  Equation  (61)  holds. 


The  Boundary  Conditions 


(a)  At  x = 


(<frx)  = (O  = (4>  ) = 0 

•*  OO  m 


(b)  At  the  wing  or  body  surface 


z W or  B 


where  9z/9x  is  the  local  slope  of  the  wing  or  body  surface 
in  the  x-direction. 

In  addition,  it  is  necessary  to  prescribe  that  the 
direct  influence  of  a disturbance  in  the  supersonic  region 
proceeds  only  in  the  downstream  direction  and  that  the 
Kutta  condition  applies  whenever  the  flow  velocity  at  the 
trailing  edge  is  subsonic  for  unique  solution. 


Integral  Equation  for  Transonic  Flow  (Non-Lifting  Case) 

Since  the  principal  object  of  the  following  analysis 
is  to  determine  the  perturbation  velocity  components  at 
any  point,  it  is  convenient  to  work  with  equivalent 
equations  for  u,  v and  w obtained  by  differentiating 
Equation  (61) : 

(i-M2)  9 + *2<fr  + d2(*>  = k i!  d2<p 
9y^  9?  9x  9x^ 
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with  respect  to  x,  y and  z,  respectively: 


( 1-M2)  ^ + ^4 


= k 


9x 


ay 


9z 


9xz 


<TT> 


(i-M2)  ^ + i^S 


9x 


3y 


9z 


k 9y9x  <T> 


(1-M2)  i-£  + L"  + a w 

ay 


9x 


= k 7 


9z 


9 z9x  V"T 


<V) 


(71) 


(72) 


(73) 


since 


9 4>  9 
9x 


2 

9x 


= u = — fu  9 
U 97  9x  C~T; 


We  will  assume  that  M<»  < 1,  i.e.  , the  free  stream  is 
a high  subsonic  velocity  though  the  local  speeds  on  the 
body  exceed  unity  in  some  areas. 

Then  normalizing  Equations  (71) , (72)  and  (73)  by 
letting 


x = x 


y = By 


z *=  Bz 


7 _ k , 

<P  - — r <P 

p 


— 9<()  k 

U = — = — n-  U 

9x  B 


— 9 <t>  k 

v = — = >j-  V 

9y  B J 


7 = d4 

9z 


k 

^w 


(74) 


where 


= \A 


m; 
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Equations  (61)  and  (71) , (72)  and  (73)  reduce  to  the 
following : 


a2(j) 


3x 


V2u  = 


2- 

V v = 


2— 

V w = 


3 2<t>  + 324>  E v2^  = 3£  3_2 <fc 
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3x  3x 


-2  -2 
3xz  * 
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<V> 


3y3x 

_3i_ 

3z3x 


-2 
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(75) 


(76) 


(77) 


(78) 


where 


,2  _ 32  . 3 2 


- + “TT  + 

3x"  3y^  3z 

is  the  Laplacian  operation  in  the  normalized  variables  x,  y 
and  z . 

Before  proceeding,  it  should  be  noted  that  the 
introduction  of  the  reduced  perturbation  velocity  component 
u permits  the  ready  recognition  of  regions  of  subsonic  and 
supersonic  velocities  and  emphasizes  the  points  at  which 
sonic  velocity  occurs.  This  becomes  apparent  from  the 
approximate  relation  between  the  local  M and  <}>  : 


1 - M2  = 1 - M2  - k 4 

CO  ~ , 


(79) 


or 


= l - 


ir 


l - M 


u = 1 - u 
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from  which  it  is  clear,  for  flows  having  subsonic  free- 
stream  M^Moo  < 1)  , that  u < 1 when  the_local  velocity  is 
subsonic,  u = 1 when  it  is  sonic,  and  u > 1 when  it  is 
supersonic . 

We  will  now  apply  Green's  theorem  to  Equation  (76). 
The  results  for  Equations  (77)  and  (78)  will  follow 
immediately.  If  a and  ft  are  any  two  functions  which, 
together  with  their  first  and  second  derivatives,  are 
finite  and  single  valued  throughout  a region  R enclosed  by 
a surface  £ , Green's  theorem,  which  follows  from  Gauss' 
divergence  theorem,  states.- 


/ 9ft  0 9(J\  _ 

(°  9K  " fi9K)dE  ‘ 


(ftV2o  - oV2ft)  dR 


(80) 


where  the  directional  derivatives  on  the  left-hand  side 
are  taken  along  the  normal  n,  drawn  inward,  to  the  surface. 


Now  in  Equation  (80)  let  ft  = u and  choose  a as  the 
fundamental  solution  l/r~  of  the  Laplace  equation  v2q  = o, 
i.e.  , J 


a 


1 

_(x-D  2+(y-n)  2+(z-c) 


m • 


Then  we  have  from  Equation  (80) 


1 9u 
9n 
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dE 


u dR 


1 

“J 


9Z 

aT3 


-2 

(^-)  dR 


since  from  Equation  (76) 
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(81) 


(82) 
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_ _ The  variable^  of  integration  in  Equation  (82)  are 

5,  n , C while  x,  y,  z are  the  coordinates  of  a point  P. 

It  must  be  observed  that  l/r3  is  singular  at  r3  = 0 and 
u is  discontinuous  at  the  shock  wave.  The  point  P and  the 
shock  must,  therefore,  be  excluded  from  the  region  R. 

A schematic  indication  of  the  body  (wing  or  fuselage)  and 
the  region  of  integration  is  shown  in  Figure  10.  The 
complete  three-dimensional  extent  of  the  body  has  not  been 
pictured.  It  suffices,  however,  to  state  that  the  surface 
£ (shown  dashed)  is  composed  of  a sphere  of  large  radius 
which  forms  the  external  boundary  of  R,  an  infinitesimal 
sphere  surrounding  P,  and  a final  surface  enveloping  the 
wing  and  body,  the  wake  (for  the  lifting  case),  and  the 
shock  waves. 

Let  us  now  apply  Equation  (82)  to  the  region  Ru 
bounded  by  the  xy-plane  and  a hemispherical  dome  of  infinite 
radius  lying  above_this  plane,  exclusive  of  the  subregions 
surrounding  P0c,y,z)  and  the  shock  waves  (Figure  10)  . Since, 
furthermore,  u may  be  assumed  to  diminish  sufficiently 
rapidly  with  distance,  the  contributions  of  the  integrals 
over  the  hemisphere  vanish  and  the  contribution  of  the 
surface  integral  over  the_small_sphere  surrounding  P may 
easily  be  shown  to  be  4ir  u(x,y,z)  when  the  radius  -*•  0 . 

We  have,  therefore,  from  Equation  (82)  the  following  result: 
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(83) 


where  th£  subscript  u denotes  conditions  on  the  upper  side 
of  the  xy-plane,  the  subscripts  1 and  2 denote  values  imme- 
diately ahead  of  and  behind  the  shock  wave,  and  S is  the 
surface  of  the  shock  wave.  The  volume  integral  is  defined 
as  follows  when  P is  ahead  of  S: 
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where  we  have  written  ip  = — — k-  (^)and  e,  , e9  are  the 

r3  z 1 z 

the  radius  of  the  small  sphere  surrounding  P(x,y,z)  and 
the  thickness  of  the  surface  surrounding  the  shock  wave 
surface  Su,  respectively.  Similarly,  when  P is  behind  S, 
the  volume  integral  is 
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If  P is  kept  fixed  in  the  upper  half  space  and  the 
region  R bounded  by  the  xy-plane  and  a hemispherical  dome 
of  infinite  radius  lying  below  this  plane  is  considered, 
it  follows  in  a similar  manner  that 
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where  the  subscript  Z denotes  conditions  on  the  lower  side 
of  the  xy-plane  and  the  volume  integral  is  defined  as 
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introducing  the  notation 


Au  = uu  - u£ 
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and  adding  Equations  (83)  and  (86) , we  have 

u(x,y,z)  = 
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where  B represents  the  wing-body  or  wing  alone  and  S repre- 
sents the  sum  of  the  surfaces  Su  and  S^;  and  R = R„  + Rj». 

The  first  integral  on  the  right-hand  side  of  Equation  (89) 
will  be  recognized  as  the  value  of  u given  by  the  linearized 
theory  of  subsonic  flow  past  wing  or  wing-body  combination. 
Denoting  this  linear  value  by  u^,  we  have 
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Equation  (89)  may  be  regarded  as  the  final  integral  equation 
for  u for  the  nonlifting  case.  But  for  numerical  computa- 
tions, it  is  advantageous  to  write  it  in  another  form  by 
integrating  the  volume  integral  twice  by  parts  with  respt 


aspect 
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to  5 , taking  proper  cognizance  of  the  definitions  given  in 
Equations  (84) , (85)  and  (37) , and  to  decompose  the  surface 
integral  over  the  shock  wave  into  components  parallel  to 
the  axes  of  the  coordinate  system.  If  then  n^ , n2  and  n3 
are  the  direction  cosines  of  the  normal  to  the  shock 
surface  drawn  inward  as  shown  in  Figure  10 , the  following 
equation  is  obtained  from  Equation  (89) : 
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where  is  given  by  Equation  (90) . 

Although  the  integral  equation  obtained  above  appears 
to  be  more  complicated  than  Equation  (89) , it  is  superior 
from  the  point  of  view  of  approximate  solutions  for  the 
following  reasons: 

(a)  The  triple  integral  of  Equation  (89)  shows  a very 
strong  influence  of  the  velocities  in  the  neighborhood  of 
P since  they  are  multiplied  by  l/r3.  This  influence  is 
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largely  nullified  in  the  triple  integral  in  Equation  (91) 
because  part  of  the  region  has  a negative  influence  and 
part  has  positive  influence.  The  predominant  influence  in 
Equation  (91)  is  furnished  by  the  term  ij2/2  standing 
outside  the  integral. 

(b)  The  contribution  of  distant  regions  is  also 
diminished  in  importance  in  the  triple  integral  of  Equation 
v^i)  since  their  influence  now  varies  as  l/r3  rather  than 
l/r3  as  in  Equation  (89). 

(c)  The  biggest  advantage  in  using  Equation  (91)  is 
that  the  value  of  the  triple  integral  in  this  equation  is 
continuous  through  a shock  wave  (due  to  the  presence  of  the 
term  u^/2  outside  the  second  derivative)  rather  than 
discontinuous  as  in  Equation  (89) . A point  of  great 
importance  in  the  approximate  solution  arises  from  the  fact 
that  the  integration  by  parts  provides  extra  terms  (those 
containing  u^/2)  in  the  integrals  along  the  shock  surface  S 
which  combine  with  those  already  present  in  such  a way  that 
the  contribution  of  these  integrals  becomes  very  small  when 
the  shocks  are  nearly  normal  waves,  as  is  usually  the  case 
for  high  subsonic  speeds. 


In  addition  to  satisfying  the  integral  equation  for  u 
given  in  Equation  (91) , the  velocity  components  on  opposite 
sides  of  shock  waves  must  be  in  accord  with  the  simplified 
relation  for  the  shock  polar  given  in  Equation  (68) , which 
can  be  rewritten  in  normalized  form,  using  Equation  (74) 
as  follows: 
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or  with  slight  manipulations  it  can  be  written  in  the 
following  two  alternative  forms: 


1 - 


ul+u2 


(u-l  - u2)2 


+ (vx  - v2)2 


+ (w ^ - w2)2  = 0 


(93) 


48 


and 
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If  the  shock  wave  is  a normal  wave  and  the  flow  is 
parallel  to  the  x-axis  (i.e.  , v]_  = v2  = wi  = w2  = 0 , 
but  ui  ^ U2> , it  can  be  seen  from  Equation  (93)  that  u 
jumps  from  1 + A immediately  ahead  of  the  shock  to  1 - A 

U1  ~ u2 

immediately  behind  the  shock,  where  A = 2 — • On  the 

other  hand,  Equation  (94)  shows  that  the  quantity  u - u^/2 
is  equal  on  the  two  sides  of  the  shock.  This  is  consistent 
with  the  fact  that  the  latter  quantity  corresponds,  in  the 
transonic  approximation,  to  the  mass  flow,  which  is 
continuous  through  a normal  shock. 

The  general  problem  of  the  three-dimensional  transonic 
flow  about  nonlifting  bodies  requires  the  solution  of 
Equation  (91)  while  taking  account  of  the  shock  relations 
given  in  Equations  (92),  (93),  and  (94).  This  is  a formi- 
dable task  well  beyond  the  reach  of  present  analysis.  We 
therefore  introduce  some  further  approximations  before 
solutions  can  be  obtained.  They  are-.  (a)  all  shock  waves 
are  assumed  to  lie  in  a plane  perpendicular  to  the  x-axis 
and  (b)  the  shock  waves  are  assumed  to  be  normal  to  the 
local  flow  direction.  The  first  assumption  corresponds  to 
setting  n2  = n^  = 0 and  the  second  leads  to  the  relations: 
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Hence  Equation  (91)  simplifies  to 
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A similar,  but  slightly  more  involved  analysis  shows 
that  from  Equations  (77)  and  (78),  the  equations  to 
determine  v and  w are 


v = vr 


w = wT 
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Before  proceeding  further,  it  should  be  observed  that 
the  solutions  of  Equations  (96) , (97)  and  (98)  must  approach 
those  of  the  linear  theory  when  <<  1,  since  u <<  1 and 
therefore  the  terms  involving  u^  °°  become  negligible, 
thereby  leaving  only 


(u,v,w)M^  <<;  1 - (ul,vl,wl)  . 

Let  us  now  write  Equation  (96)  as  follows: 
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where 
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Although  I is  a function  of  u and  is  therefore  unknown, 
it  is  informative  to_rewrite  Equation  (100)  by  solving  for 
u in  terms  of  I and  u. , thus 


u = 1 ± /I  - (2ul  - 1)  = 1 ± /I  L 


(102) 
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where 


L = 2ul  - 1 


from 

u 


(103) 


or 

supersonic.  A change  in  sign  at  a point  where  the  radical 
is  zero  corresponds  to  a discontinuous  jump  in  velocity. 

As  pointed  out  following  Equation  (94) , such  discontinuities 
correspond  to  normal  shock  waves  and  are  permissible  when 
they  proceed  from  supersonic  to  subsonic  velocities  or  from 
plus  to  minus  sign  in  Equation  (102) , when  progressing  in 
the  flow  direction.  Discontinuities  in  the  reverse 
direction  are  inadmissable  since  they  correspond  to 
expansion  shocks,  a phenomenon  which  violates  the  second 
law  of  thermodynamics. 

The  values  of  u^,  and  hence  L,  can  be  calculated  for 
any  given  body  and  are  generally  characterized  by  certain 
regions  in  which  ul  > 0_and  other  regions  in  which  ul  < 0. 
The  absolute  values  of  u increase  continuously  with 
increasing  Mro  and  the  maximum  positive  values  may 
considerably  exceed  unity  as  sonic  velocity  is  approached 
in  the  free  stream.  Not  very  much  can  be  stated  at  this 
point  about  the  values  of  I,  except  that  they  depend  on  the 
distribution  as  well  as  magnitude  of  u and  that  the  above 
inequality  must  be  satisfied. 


Several  interesting  observations  can  be  made 
Equation  (102).  First  of  all,  for  real  values  of 

I > L . 


Furthermore,  the  choice  of  the  plus  or  minus  sign 
determines  whether  the  local  velocity  is  subsonic 


The  Two-Dimensional  Flow  Past  Circular  Arc  Airfoil: 

The  Oswatitsch-Spreiter-Alksne  Solution  of  the  Integral 
Equation 

The  starting  point  for  the  solution  of  the  full 
integral  equation  derived  for  the  three-dimensional  cases, 
would  be  to  examine  the  methods  needed  for  two-dimensional 
flows  given  by  a number  of  authors:  Oswatitsch,  Gullstrand, 
Spreiter  and  Alksne , Nixon,  and  Norstrud. 


The  analysis  of  the  two-dimensional  case  is  based  on 
the  following  equation: 


—2  /—  - 


—2  /-=-  -s 


'77  2 /—  — \ 2 . 


u(x,z)  =uL(x,z)  . ||  u (|,0  Lfr’O  r(g-0 


[ (x-o2+(z-o2] 


dCdc 


Approximate  solutions  of  this  equation  could  conceivably 
be  worked  out  numerically  by  starting  with  a two-dimensional 
grid  of  suitably  selected  values  of  u and  iterating  until_ 
convergence  is  obtained.  If  the  first  approximation  for  u 
is  taken  to  be  the  results  given  by  incompressible  or 
linearized  subsonic  flow  theory,  convergence  will  be  obtained 
only  when  Mra  is  sufficiently  small  and  the  flow  is  subsonic 
everywhere.  Oswatitsch  suggested  that  mixed  flow  fields 
containing  shock  waves  can  be  obtained  if  the  starting 
u-distribution  contains  shock  waves.  The  idea  then  is  to 
start  with  a reasonable  guess  for  u,  being  sure  to  include 
a proper  discontinuity  complying  with  the  shock  relations 
of  Equation  (95),  then  proceed  to  the  solution.  It  is  not 
necessary  to  be  highly  accurate  in  the  initial  guess  for  u. 

A further  simplifying  assumption  is  made  by  Oswatitsch 
and  others  regarding  the  variation  of  u in  the  z-direction, 
which  has  the  effect  of  reducing  the  double  integral  in 
Equation  (104)  to  a single  integral.  Oswatitsch  suggested 
the  following  relation: 


u(x,z)  = 


Ug (x , 0) 


u + h 


(105) 


where  ug(x,0)  is  the  value  of  u on  the  body  as  taken  in  the 
linearized  theory  and  b is  a function  of  x so  chosen  that 

the  irrotationality  condition  ^ is  fulfilled  at  z = 0. 
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SECTION  V 


THE  INTEGRAL  EQUATION  FOR  THE  THREE-DIMENSIONAL 
LIFTING  WING- BODY  COMBINATION 


We  start  with  the  nonlinear  partial  differential 
Equation  (61) : 


(1  - M ) <p  + <p  + (p  — k 0 <J> 
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where 
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The  change  of  type  of  the  equation  is  an  essential 
feature  of  transonic  flow  and  the  regions  where  this  occurs 
are  recognized  by  the  sign  of  the  coefficient  of 
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Equation  (106)  is,  of  course,  valid  only  in  regions 
where  the  necessary  derivatives  exist  and  are  continuous. 
Across  a shock  surface  the  normal  component  of  velocity  is 
discontinuous  while  the  tangential  component,  and  therefore 
<J) , is  continuous.  The  necessary  relation  follows  from  the 
classical  expression  for  the  shock  polar  given  in  Section  IV. 
This  approximates,  in  small  disturbance  transonic  theory, 
to  [see  Equation  (68)]: 
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where  the  subscripts  1 and  2 refer  to  conditions  ahead 
of  and  behind  the  shock. 


Equations  (106)  and  (109)  are  usually  applicable  to 
the  case  where  the  x-axis  is  parallel  to  the  free-stream 
at  infinity  and  are  good  approximations  where  the  coordinate 
system  is  rotated  slightly.  Choosing  the  body-axis  system, 
whrre  x-axis  is  aligned  in  the  longitudinal  axis  of  the 
wing  and  body,  the  relation  between  the  total  velocity 
potential  4>(x,y,z)  and  the  perturbation  velocity  potential 
4>(x,y,z)  is  approximated  by: 


4>(x,y,z)  = Uoo(x  + az)  + <j>(x,y,z)  (110) 


where  a is  the  angle  of  attack. 

The  expression  for  the  pressure  coefficient  Cp  is  not 
invariant  with  respect  to  small  rotations  of  the  coordinate 
system.  In  body  axes,  the  proper  expression  is 


Cp  - - U:  (fx  + o+z)  - jfl  <*y2  + »z2>  (U» 


The  Boundary  Conditions 

The  condition  at  infinity  yields 


<t>(°°)  = Uoo(x  + az) 


or  that 


<K“)  = 0 . (112) 


An  exception  to  this  statement  occurs  in  the  vicinity  of 
the  wake  at  great  distances  behind  the  wing,  but  no 
complications  arise  due  to  the  relative  smallness  of  this 
region.  The  conditions  at  the  airplane  surface  result  in 
the  relation 
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where  n]_,  n2  and  W3  are  the  direction  cosines  of  a normal 
to  the  airplane  surface  with  respect  to  the  x,  y,  and  z 
axes,  respectively.  This  relation  is  too  general.  For 

our  needs  n^  <<  1 or  (n^+no)  . i.e.,  slender  bodies  and 

wings  are  considered.  Then  "'Equation  (113)  reduces  to 
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where  n is  the  normal  to  the  curve  bounding  a cross  section 
in  a plane  normal  to  the  x-axis. 

We  will  now  apply  Green's  theorem  in  two  different 
forms,  one  for  S 1,  and  the  other  for  M > 1.  Let 
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Case  1.  Subsonic  Freestream  Flow:  5 1 

Equation  (106)  can  be  rewritten  as 

8 2 

|2*«  + *yy  + *zz  ’ k*x*xx  = k K <-J->  U16) 


and  the  corresponding  expression  of  Green's  theorem  is 
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where  ft  and  4>  are  arbitrary  functions  and  L is  the 
Prandtl-Glauert  operator: 
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and  is  the  derivation  along  the  conormal  defined  by 
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where  ni,  n2 , n3  are  the  direction  cosines  of  the  normal 
to  the  surface  drawn  into  the  region  V. 

Now  choose  ip  = , where 


- 1/2 
0=  (x  - O2  + 62((y  - n)2  + (Z  - 02f|  (120) 


so  that  L(ip)  = 0 and  let  ft  = <|>  , the  perturbation  velocity 
potential,  in  Equation  (117).  Then  from  Equations  (116) 
and  (117)  we  obtain 
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In  these  equations  the  running  coordinates  in  the 
integrations  are  n , i,  and  <p  is  to  be  calculated  at 
point  P (x , y , z) . 

If  Equation  (121)  is  applied  to  the  infinite  region 
V consisting  of  a sphere  of  large  radius  which  forms  the 
external  boundary  V,  a sphere  of  infinitesimal  radius 
surrounding  the  field  point  P(x,y,z),  and  a final  surface 
enveloping  the  wing-body,  its  wake,  and  its  shock  waves, 
we  obtain  the  following  expression: 
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where  B+W  denotes  the  wing-body  surface  and  its  wake  and 
^1  fnd  ^2  are  the  shocks,  upstream  and  downstream  shock 
surfaces,  respectively.  On  the  shock  surface  A the 
conormals  are  directly  opposed  on  the  upstream  and  down- 
stream faces  A]_  and  A2.  On  the  body  itself  and  the  wake, 
as  mentioned  earlier,  we  may  approximate 
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If  now  the  triple  integral  in  Equation  (122)  is 
integrated  by  parts  twice  with  respect  to  x,  the  resultant 
form  of  Equation  (122)  is 
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Equation  (123)  is  of  particular  interest  because  the 
integrals  over  the  shock  surfaces  may  be  shown  to  vanish  as 
Ln  no^lifting  case.  Finally,  if  we  neglect  the  term 
n]#^  in  comparison  with  3<j>/3v  , the  simplified  integral 
equation  is  similar  to  the  nonlifting  case: 


57 


£!3 


'(x,y  ,z)  = -^ 


riM.  _ 

'o  3n 

B+W 


3 1^rlV 
3Ha)di: 


k 

5¥ 


f[ 

V 


1 

7 


<^dV 


(124) 


The  first  integral  on  the  right-hand  side  of 
Equation  (124)  is  the  expression  for  <J>(x,y,z)  in 
linearized  theory  and  the  spatial  integral  is  a 
contribution  brought  about  by  the  nonlinear  term  of  the 
basic  differential  Equation  (106).  It  is  of  interest 
to  remark  that  a straightforward  derivation  ignoring 
the  existence  of  the  shock  waves  leads  to  the  same 
equation  as  Equation  (124).  For  the  majority  of  cases 
of  practical  interest,  it  can  be  shown  that  compensating 
terms  arise  and  the  discontinuity  surfaces  are  taken 
care  of  by  a formal  development  that  ignores  the 
existence  of  these  surfaces.  However,  if  it  is  not 
possible  to  ignore  so  completely  the  existence  of  the 
discontinuity  surface,  Equation  (122)  will  be  preferable 
to  Equation  (124). 

Since 
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Equation  (124)  may  be  written  in  the  final  forms  as 
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Introducing  the  notation 
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and  the  variable 
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and  remembering  that  <(>  is  continuous  at  the  shock 
surface  the  more  exact  Equation  (122)  may  be  written 
in  the  following  two  alternative  forms : 


4>(x,y  ,z)  = <J>L(x,y  ,z) 
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or 


(Kx.y.z)  = <J>L(x,y,z)  - ^ ||A  [|^]  WdE  - ^ ^ |f 


> Vl 
raT  -i-J wdv 


(128) 


The  perturbation  velocity  components,  u,  v and  w, 
are  obtained  by  differentiating  cj>  with  respect  to  x,y,z, 
respectively  from  either  Equation  (125/126)  or 
Equation  (127/128). 

Consider  Equation  (125)  for  the  x-derivative . After 
isolating  the  singularity  at  the  field  point  by  intro- 
ducing the  limits  5 = x ± e and  using  the  rule  of 
differentiating  under  the  sign  of  integration  with 
variable  limits : 
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In  the  limit  as  e-*0  the  influence  function  in  the 
integrand  of  the  double  integral  is  a two-dimensional 
Dirac  delta  function  at  the  point  n = y , C = z and  of 
strength  2 tt / g2  . The  expression  for  u then  becomes 
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u(x,y,s)  =uL(x,y,s)  +^.-ai%y,z>  . k.  |J|  £(i^l 
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(129) 
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where  we  have  made  use  of  the  fact  that 

_9_  1 _ _3_  1 

ay  a an  a 
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and  the  double  integral  tends  to  zero  in  the  limit  for 
all  values  of  E,  and  t,  . 


Similarly  it  is  easily  seen,  by  differentiating 
Equation  (125)  with  respect  to  z,  that 


w(x,y,z)  = wL(x,y,z)  - ^ 
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Case  2.  Integral  Equation  for^ M 1 1 


If  the  undisturbed  flow  at  infinity  is  supersonic, 
i.e.,  Mro  > 1,  Equapion  (106)  can  be  written  as 


6$  + 4>  + 6 = k d>  d>  = k ^ x 

xx  yyy  Tzz  vxvxx  3x 


(132) 


and  the  corresponding  expression  of  Green's  theorem  is 
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where 
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and 
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The  derivation  of  the  integral  equation 
in  Case  1 would  suggest  that  ip  be  replaced  by 

[(x-O2  - 32(y-n)2  - ■'•^2  but  this  leads  to  the 

introduction  of  a finite-part  technique  in  the  integration. 
For  the  initial  stages  of  the  analysis,  ip  will  be 
identified  with  the  fundamental  solution 


cos 


X - ? 

- ~ TO 

B[(y-n)z  + (z-Oz] 


of  L(i p)  = 0 used  by  Volterra  (see  Reference  14,  page  ). 
Then  from  Equations  (132)  and  (133)  , we  have 


♦ ^ dr 

8 V* 


|j|  o»L  (4>)dV 
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(136) 


We_have  now  to  choose  properly_the  three-dimensional 
region  V and  its  enclosing  surface  I . Discontinuities 
in  the  velocity  components  are  again  to  be  taken  into 
consideration  at  the  shock  waves.  Furthermore,  the 
fundamental  solution  of  Volterra  becomes  infinite  at 
n = y , C = z , that  is,  everywhere  along  the  line  passing 
through  the  field  point  P and  parallel  to  the  x-axis. 
Figure  11  indicates  the  disturbance  field  of  the  body 
as  well  as  V and  Z (dashed  lines) . The  bow  shock_fixes 
the  foremost  extent  of  the  disturbance  field  and  Z lies 
adjacent  to  it  and  other  possible  shock  surfaces  as  well 
as  the  surface  of_the  body  and  its  wake.  The  downstream 
limits  of  region  V are  fixed  by  the  forecone  with  vortex  P 
and  determined  explicitly  by  the  relation 


2 2 1/2 

x - K = e[(y-n)  + (z-Oz] 


(137) 


Bow  Shock 


Figure  11.  Region  of  Integration  for  £ 1. 


The  inner  boundary  of  V is  the  cylindrical  surface 
of  infinitesimal  radius  given  by  the  relation 


(y-n)2  + (z-O2  = e2 


The  conormal  derivative  is  defined  by  Equation  (135) . 
On  the  infinitesimal  cylinder  its  direction  is  parallel 
to  that  of  the  normal  to  the  surface , and  on  the  forecone 
from  P the  conormal  is  directed  along  the  surface  itself. 
We,  therefore,  obtain  from  Equation  (136)  the  following 
expression  for  <t>(x,y  ,z)  : 
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where  integrals  over  the  surface  of  the  body  and  wake  are 
denoted  by  t , over  the  two  sides  of  the  shock  surfaces  by 
A]_  and  \2 , and  over  the  enclosed  volume  by  V.  In  each 
case,  only  that  portion  of  the  surface  or  volume  lying 
within  the  forecone  of  P is  included  in  the  integrals. 

The  surface  integrals  over  the  forecone  itself  vanish 
because  w and  9w/ 9v  are  zero  on  it  (Reference  14) . 

Integration  by  parts,  in  the  volume  integral  in 
Equation  (138),  leads  to  the  relation: 
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Equation  (139)  is  the  form  for  H,,  _ 1 analogous  to 
Equation  (123)  for  < 1 and  on  the  body  and  wake  surface 
involves  the  approximation 


„ 3 1 3 
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As  before,  it  can  be  shown  from  the  shock -wave 
relation  that  the  combined  integrals  over  the  surfaces 
and  A2  vanish.  Therefore 
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Since 
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Equation  (140)  may  be  written  in  the  following  two 
alternative  forms: 


<J>(x,y,z)  = <{>L(x,y  ,z) 
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or 
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<J)(x,y  ,z)  = (J»L(x,y  ,z) 
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Similarly,  Equation  (138)  may  be  given  in  the 
following  two  alternative  forms  using  the  notation 
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(145) 


where  <j>^  is  the  value  from  the  linearized  supersonic  theory. 


Conditions  at  the  Shock  Surface 

It  is  of  some  interest  to  study  Equation  (125)  or 
Equation  (126)  as  the  field  point  approaches  a discontinuity 
surface  and  to  discover  the  mechanism  by  means  of  which 
these  basic  equations  furnish  the  velocity  jumps  associated 
with  the  shock  waves  in  the  field.  The  surface  of  the  wave 
can  be  replaced  locally  by  a planar  element  and  a new 
coordinate  system  X,Y,Z  introduced  with  the  origin  fixed 
at  the  intersection  of  the  line  n = y , t,  = z and  the  shock 
surface.  Then  the  point  P,  at  which  conditions  are  to  be 
calculated  close  to  the  surface  of  the  wave . has  the 
coordinates  (X,0,0)  and  the  planar  surface  of  the  wave  is 
given  by  the  linear  relation 


aX1  + bYx  + cZL  = 0 . 


(146) 


67 


We  assume  that  u in  Equation  (125)  is  composed  of  a 
continuous  part  and  a discontinuous  part  that  has  the 
constant  value  u?  ahead  of  and  u^  behind  the  shock. 
Equation  (125)  then  yields 

Urn  [u(X  ,0,0)  - u ( X , , 0 , 0 ) ] 

X->0 


„ . k , 2 2n  9 

■ tlm  (uru2>  3x 


a dYx  dZ1 

[ (aX+bYj+cZ^2  + a2B2Y12+a2B2Z12] 


(147) 

where  the  double  integral  extends  over  the  region  of  the 
discontinuity.  If  the  differentiation  with  respect  to  X 
is  now  carried  out  within  the  integral  signs  and  X allowed 
to  approach  zero,  we  will  obtain 


u-.  - u9  = £im  [u(X  ,0,0)  - u(X,  ,0,0)] 
L Z X-0 


2 2 
ul-u2 


ag  +b  +c" 


(148) 


It  can  easily  be  shown  that  Equation  (148)  agrees  with 
the  result  given  by  Equation  (68)  , the  shock  polar  condi- 
tion , for  M £ 1 . 

' no 


SECTION  VI 


SOLUTION  OF  THE  TRANSONIC  THREE-DIMENSIONAL 
INTEGRAL  EQUATION  FOR  THE  NONLIFTING  CASE 


We  have  described  the  difficulties  and  complications 
involved  in  the  solution  of  two-dimensional  transonic 
nonlinear  equation.  As  a consequence  of  the  obstacles 
encountered  in  the  two-dimensional  case,  attempts  at  the 
solution  of  the  nonlinear  equation  with  three  space 
variables  are  very  few.  In  fact,  only  two  references 
can  be  cited  which  can  be  regarded  as  a proper  three- 
dimensional  analysis.  The  first  one  is  the  work  by  Alksne 
and  Spreiter  (Reference  15)  which  is  actually  an  extension 
of  the  local  linearization  method  previously  developed  for 
planar  (Reference  16)  and  axisymmetric  flow  (Reference  17) . 
This  has  been  summarized  below.  The  second  work  is  by 
Norstrud  who,  for  the  first  time,  attempts  to  solve  the 
nonlinear  integral  equation  by  analytical-cum-numerical 
technique.  This  needs  detailed  review  and  discussion. 


1.  THE  ALKSNE -SPREITER  METHOD  OF  SOLUTION  OF  THE  TRANSONIC 
THREE-DIMENSIONAL  EQUATION  FOR  NONLIFTING  WINGS 

Essentially  this  method  is  an  extension  of  the  local 
linearization  technique  so  successfully  employed  by  the 
authors  for  two-dimensional  thin  airfoils  in  transonic 
flows  (Reference  16) . 

The  axis  system,  as  shown  in  Figure  12,  consists  of 
the  Cartesian  coordinates  with  x-axis  parallel  to  the  free- 
stream  and  with  the  origin  at  the  west  forward  point  of  the 
wing.  The  wing  sections  are  given  by  z = Z(x,y)  and  the 
wing  planform  is  described  by  y = -S^(x)  and  y = S£(x)  as 
shown.  For  all  practical  wings  we  will  take  |S9  | = 

I S lmax I = b/2>  the  semispan.  max 

If  <|>  is  the  perturbation  velocity  potential,  whose 
gradient  yields  the  perturbation  velocity  components  u,  v 
and  w,  parallel  to  the  x,  y and  z axes,  the  problem  of 
transonic  flow  around  thin  finite  wings  can  be  studied  by 
using  the  usual  nonlinear  equation. 


( 1 - M2  ) 0 + 
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(149) 
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The  appropriate  boundary  conditions  are: 

(1)  u,v,w  -*•  0 far  ahead  of  the  wing,  and 

(2)  the  flow  be  tangential  to  the  wing  surface. 

The  latter  condition  is  satisfied  if 

(*  ) - u 3Z(*'y)  ■ (150) 

z z-0  " 9x 


In  addition  to  the  above  condition  it  is  necessary 
to  take  proper  account  of  the  difference  in  regions  of 
influence  and  dependence  in  the  subsonic  and  supersonic 
portion  of  the  flow  field.  In  general,  additional 
equations  are  needed  for  flows  across  the  shock.  Alksne 
and  Spreiter  did  not  consider  them  for  the  following 
reasons : 

(1)  In  the  first  place,  attention  was  restricted  to 
Meo  very  close  to  1 so  that  the  flow  in  the  vicinity  of  the 
wing  is  essentially  the  same , by  virtue  of  the  Mach  number 
freeze , as  at  = 1 . 

(2)  Secondly,  they  assumed  that  the  shock  stands  at 
the  rear  of  the  wing  where  its  influence  cannot  extend 
onto  the  wing  surface. 


Approximate  Solution 

First,  the  coefficient  of  <f>x  is  replaced  by  the 
symbol  X as  follows: 


X = 


M^(y+1) 


U 


xx 


5 k i > 0 


so  that  Equation  (149)  becomes: 


t,  +<|)  -X*  = -(l-  M ) <(> 

yy  vzz  yx  v ®/Yxx 


(151) 


(152) 


They  have  restricted  their  attention  to  only 
accelerated  flows  for  which  X > 0.  In  the  initial  stages 
of  the  analysis  it  is  assumed  that  X is  a constant. 
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The  resulting  Equation  (152)  is  linear  for  X = constant 
and  is  oi  elliptic,  hyperbolic  or  parabolic  type  depending 
on  whether  1.  Use  of  the  boundary  conditions  given 
above  anc  the  application  of  Green's  theorem  to  Equation 
(149)  leads  to  the  following  equation  for  u if  shock 
conditions  are  not  taken  into  account: 


up(x,y,z)  = - 
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(153) 


where 


f«  = - 


At  = 1 , f p = 0 and  up  can  be  calculated  for  any 
given  value  of  X for  a known  profile  and  planform.  At 
other  Moo  Equation  (153)  is  an  integral  equation  whose 
solution  remains  to  be  found.  For  Moo  near  1,  however,  fp 
is  very  small  and  <t>frp  can  safely  be  replaced  by  X/k  from 
Equation  (151).  The’triple  integral  can  then  be  evaluated 
giving  us  the  following  result : 
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dn 


Next,  the  value  of  k<J>xx  at  the  point  x,y,z  is  restored 
in  place  of  X in  Equation  (154) . Then  for  each  value  of 
y and  z the  result  is  a first-order,  nonlinear,  ordinary 
differential  equation  for  u as  a function  of  x.  Dropping 
the  subscript  P and  writing  k<t>xx  as  k u'  , we  obtain  from 
Equation  (154): 
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It  is  convenient  to  express  Equation  (155)  in  terms 
of  reduced  variables  indicated  by  the  transonic  similarity 
rule  and  defined  as  follows: 
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Equation  (155)  thus  becomes 


u(x,y,z)  = 
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In  cases  where  interest  is  confined  to  conditions 
on  the  wing  surface,  u can  be  evaluated  with  z = 0 and 
Equation  (157)  can  be  simplified  to  become 


u(x,y,0)  = 


u' =const 


x S2(Z) 
dC 
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3(Z/t)/3g 

x-Z 
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(158) 


Properties  of  Transonic  Flow 

From  Equation  (158)  it  can  be  seen  that  for  affinely 
related  planforms  (fixeds(x^)  and  affinely  related 
profiles  (fixed  3(Z/t)/9x),  u at  a given  x and  y varies 
only  with 


2 

M -1 


= 


[M^,(Y+1)t] 


T/J  • 


Thus  the  present  results  follow  the  transonic 
similarity  rule.  Note  also  that  s has  the  same  parametric 
form  as  the  aspect  ratio  parameter  customarily  used  for 

the  transonic  similarity  rules,  i.e.,  A = m£(y+1)t]^^  A, 

where  A is  the  aspect  ratio. 

Furthermore ,_u'  in  the  exponent  in  Equation  (158)  can 

be  replaced  by  d(u^°°)  for  a given  value  of  y and  it 

dx 

follows  for_wings_of  affinely  related  geometry  that  u+5<» 
at  a given  x and  y does  not  vary  with  Thus  the  results 

given  by  Alksne  and  Spreiter  contain  the  Mach  number  freeze. 
Notice  also  that  there  is  a coupling_of  aspect  ratio, 
thickness  ratio,  and  Mach  number  in  A,  such  that  in  order 
to  maintain  a fixed  A for  a wing  of  given  thickness  ratio 
as  M»  varies  it  is  necessary  to  vary  the  aspect  ratio 
likewise.  The  variation  of  aspect  ratio  with  M^  is  not 
great  in  the  range  for  which  the  freeze  can  be  expected  to 
apply. 


74 


A 


For  each  value  of  y and  z Equation  (157)  is  a first 
order  ordinary  differential  equation  and  its  solution 
requires  the  specification  of  an  initial  or  a boundary 
condition.  The  simplest  case  is  one  for  which  a value  of 
u is  known  at  some  value  of  x for  every  spanwise  station 
and  z = 0.  This  is  the  case  for  a wing  of  wedge  profile, 
for  which  the  sonic  velocity  occurs  at  the  shoulder. 

In  general,  however,  there  is  no  point  at  which  u is  known 
a priori  and  some  alternative  condition  must  be  imposed 
in  order  to  determine  a unique  solution.  It  happens  in 
all  cases  considered,  that  is,  two-dimensional  flow  around 
an  airfoil,  axisymmetric  flow  around  a body  of  revolution, 
and  now  three-dimensional  flow  past  a wing,  that  the  equation 
contains  a singular  point  through  which  pass  an  infinity  of 
integral  curves.  In  each  case  only  one  is  analytic.  Thus 
the  requirement  of  analyticity  is  sufficient  to  determine 
a unique  solution. 

Mathematically,  this  leads  to  the  following  condition: 
Writing  Equation  (158)  as 

U = F(x,U' ;y)  (159) 

where 

U = u + £ 

^00 

and 

d(u+£  ) 

U*  = u'  = 

dx 

the  equation  will  be  investigated  for  a fixed  value  of  y. 

From  Equation  (159) , we  have 


U'  = El  + u"  £l  (160) 

3x  3u' 

where,  in  general,  3F/3x  and  3F/3U'  are  functions  of  x 
and  U'  . 

At  a point  where  3F/3U'=  0 we  can  find  U'  provided 
that  U"  is  finite,  but  this  is  assured  if  the  solution  is 
analytic.  _With  U'  known  at  a point,  it  is  possible  to 
calculate  U at  that  point  using  Equation  (159) . This  is 
not  enough  to  determine  a unique  solution  because  U"  can 
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have  any  finite  value.  However,  additional  derivatives 
can  be  taken  and  in  each  case  the  unknown,  but  finite, 
higher  derivative  of  U will  be  multiplied  by  9F/9U' 
which  is  equated  to  zero.  Thus  all  derivations  can  be 
determined  in  principle  at  the  singular  point. 

The  location  of  the  singular  point  and  the  velocity 
g^idient  at  that  point  are  then  obtained  by  solving 
simultaneously  the  two  equations: 


5F 
9U ' 


0 and  U'=  ^ 
9x 


Special  Cases 

Alksne  and  Spreiter  have  considered  several  special 
cases  and  deduced  closed  form  expressions  for  Equation 
(157)  or  Equation  (158). 


(1)  Central  Region  of  High  Aspect  Ratio  Wing 

By  letting  S^(x)  = ^(x)  = 00  and  considering  Z a 
function  of  x only  across  the  entire  span,  the  following 
asymptotic  form  of  Equation  (158)  results: 
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If  the  sonic  point,  x*,  where  u + = 0,  can  be 

found  or  is  known,  the  general  solution  may  be  written  as 
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For  a smooth  profile  for  which  u + ^ is  not  known 
a priori  at  any  point,  the  singular  point°°is  found  to  be 
the  value  of  x for  whicu 
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Using  this  value  of  x*  in  Equation  (162)  a unique 
solution  is  determined. 

Theoretical  pressure  distributions  computed  for 
several  airfoils  using  the  above  equations  are  in  good 
agreement  with  experimental  pressure  distributions  at 
M = 1 (Reference  16) . 


(2)  Tip  of  High  Aspect-Ratio  Wing 

_P lacing  the  origin  at  the  tip  so  that  y = si  = 0 
and  S2(x)  can  be  set  equal  to  infinity,  Z can  then  be 
considered  a function  of  x alone  across  the  entire  span. 
Equation  (158)  then  becomes 


u + £ - - 

' rr\ 


X 


TT  U 


_d_ 

dx 


d(Z/t)/d£ 
~ Z — ZTT77 
(x  - 0±U 


(164) 


This  is  the  same  result  as  Equation  (161)  with  a 
factor  1/2  in  Equation  (164)  which  cannot  affect  the 
position  of  the  singular  point.  It  follows  therefore 


(3)  Small  Aspect-Ratio  Wing 

In  Equation  (155)  let  the  integration  be  performed  in 
two  steps-.  0 to  x - e and  then  from  x - e to  x,  where  e 
is  a small  quantity. 


x-b  *: 

iu,y,z)_  1 3 1 ^ f 

~U 2¥  3xJ  , J dM 

00  J u =const  _JC 


.x-e  So (5) 


9Z/3£ 


-ku'[(y-n)2+z2] 


■S1(€) 


S,(0 


X 0 9 

I «j  3 


x-e  -S-j^  (C) 


1 - M * 

00 

M 2(y+l) 


(166) 


Since  the  point  £ = x is  excluded  from  the  first  term 
on  the  right-hand  side,  the  exponential  term  for  a small 
aspect  ratio  wing  can  be  replaced  by  1.  In  the  second 
integral  we  replace  3Z/3£  = 3Z/3x  in  the  small  interval 
(x-e,x)  and  obtain  the  simpler  equation: 

' x-e  S2(£) 

u _ 1 3 f d£  \ 3Z  , 

U 27  3x  , ‘ J 5TT£  3£  dn 

J u =const  ■'q  -§l(£) 


S2(x) 
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where  S(x)  is  the  cross-sectional  area  of  the  wing 
perpendicular  to  the  x-direction  if  Z = 0 at  the  leading 
and  trailing  edges  at  any  spanwise  station  and  the 
?- integration  of  the  second  term  on  the  right  can  be 
performed  approximately  by  retaining  only  the  first 
logarithmic  term. 

The  desired  asymptotic  form  for  a wing  of  vanishing 
aspect-ratio  becomes 
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where  the  Euler's  constant  C = 0.5772156 


Remarks  on  the  Alksne -Spreiter  Method 

The  method  is  elegant  and  simple  to  use  with  fair 
degree  of  accuracy  but  it  is  very  limited  in  scope  as 
regards  the  flow  conditions  and  the  wing  geometry.  The 
following  are  some  of  the  restrictions : 

(a)  The  wings  must  have  sharp  leading  edges  with 
zero  lift.  Randall  (Reference  18)  has  given  an  extension 
applicable,  under  certain  circumstances,  to  round  nosed 
two-dimensional  airfoils  at  angles  of  attack. 
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(b)  The  Mach  number  should  be  very  near  one  and  the 
wing  planform  and  profile  should  be  such  that  at  all  points 
the  flow  accelerates  along  the  chord. 
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(170) 


The  physical  interpretation  of  Equation  (169)  is  that 
of  a mutual  interference  between  velocity  points  in  a 
perturbed  flow-field  u(x,y,z)  where  the  linearized  field 
uL(x,y,z)  designates  the  datum  for  zero  perturbation.  The 
kernel  of  the  triple  integral  serves  as  the  influence 
measure  and  it  can  be  seen  from  Equation  (169)  that  this 
function  is  largely  dependent  on  the  double  derivative  in 
the  x-direction  of  the  inverse  of  the  distance  r3  between 
influencing  points.  This  means  that  the  strongest  influence 
at  a point  P(x,y,z)  in  the  flow  field  will  come  from  field 
points  which  are  located  directly  upstream  or  downstream  of 
P.  The  spanwise  effects,  however,  are  of  major  concern  in 
describing  spatial  transonic  flow.  Here  the  compressibility 
of  the  flowing  gas  tends  to_freeze  the  available  stream  tube 
area  and  the  density  flux  pu  must  be  balanced  with  neighbor- 
ing stream  tubes.  Norstrud  assumes  the  following  exponential 
functional  relationship  to  extract  the  spanwise  influence 
from  Equation  (169)  : 
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u(x,y,z)  = u(x,y ,0)e“Z^r  (171) 

which  gives  a relation  between  the  velocity  in  the 
xy-plane  of  the  reduced  wing  and  the  velocities  directly 
above.  The  realation  [Equation_(171) ] satisfies  correctly 
the  boundary  condition  of  zero  u-disturbance  at  infinity 
(z  -*■  00)  . The  parameter  r = r(x,y)is  defined  using  the 
irrotationality  condition  and  the  linearized  value  uj/x.y.O) 
over  the  wing  which  automatically  satisfies  the  tangential 
flow  condition: 


r(x,y)  = - 
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1 x ,y , 0 

Using  the  irrotationality  condition 
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we  have  from  the  linearized  boundary  condition 


uT (x,y ,0) 

r(x,y)= 
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where  Z represents  the  reduced  ordinates  of  the  wing  section 
at  station  y. 

Using  the  tables  of  the  Laplace  transform,  the 
Equation  (169)  can  be  rewritten,  for  x - £ = 0,  as 


— i_2  — — — — ___ 

u(x,y,z)  = j-u  (x,y,z)  +uL(x,y,z)  -e 


,Z  -H» 


u2(£,n,0)E3(|y-ri|  ;r)d£dn 


(174) 


The  spanwise  influence  function  E^  = E3 ( | y -n | ; r)  is 
found,  for  the  argument  a = 2|y-n|/r  (from  Reference  19) 
to  be 
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For  small  values  of  the  argument,  we  have 


SiX 


X X3  x X5 

TTiy  ' TTTT  + 5T5T  ~ • 


C = Euler's  constant  = 0.57721  .... 

and  hence  yields  a logarithmic  singularity  for  |X|  -*•  0 
The  rapid  decrease  in  spanwise  influence  is  typical  in 
spatial  flows.  The  lateral  influence  function  Eo(o)  is 
shown  plotted  in  Figure  13.  J 


Figure  13.  Lateral  Influence  Function  E~(o). 


Norstrud's  Approximate  Method  of  Solution  of  Equation  (169) 


The  symmetrical  wing  of  arbitrary  planform  is  repre- 
sented by  a number  of  rectangular  panels  in  the  xy-plane. 

The  distribution  used  for  these  panels  is  shown  in  Figure  14. 


>*»  tits 


Figure  14.  Approximate  Representation  of  the  Wing. 


The  fundamental  point  of  view  adopted  by  Norstrud  in 
the  definition  of  the  influence  region  for  a projected 
point  of  interest  (x^y^.O)  on  the  wing  for  the  evaluation 
of  the  double  integral  in  Equation  (174)  is  to  collect  the 
effects  of  all  field  points  located  at  a common  streamwise 
station  (x  = x ,x,,  . . . , x ) and  assume  this  total 
influence  to  be1representednat  the  point  (x,y0,o)  by  a 
mean  value  obtained  from  two-dimensional  analysis.  These 
influences  which  may  constitute  the  solution  from  a strip 
method  are  supplemented  by  the  effects  of  laterally  located 
points,  i.e.  , for  x - xQ  and  ytyp.  The  influence  pattern 
is  illustrated  in  Figure  14  by  the  shaded  area  and  the 
resulting  influence  matrix.  [I]  has  the  general  structure 
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The  influence  coefficients  ev"  ’ . are  integral  functions 

1 » J 

of  the  defined  influence  functions  E2  and  E3  for  the  case 

where  £ = j and  k = i(£  ^ j) , respectively.  If  £ # j and 
1-  0 

k + i,  e. = 0.  Equation  (174)  can  now  be  written  in  the 
^ 1 J 

matrix  form 
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(177) 


and  the  problem  has  been  reduced  to  that  of  solving  a system 
of  nonlinear  algebraic  equations.  Norstrud  used  the  Newton- 
Raphson  method  of  successive  approximations  to  solve  the 
system  [Equation  (177)].  This  utilizes  a linear  system  at 
each  iteration  step. 
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